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Abstract 

We use recently calculated next-to-next-to-leading (NNLO) anomalous di- 
mension coefficients for the n = 1, 3, 5, . . . , 13 moments of the xF^ structure 
function in vN scattering, together with the corresponding three-loop Wil- 
son coefficients, to obtain improved QCD predictions for the moments. The 
Complete Renormalization Group Improvement (CORGI) approach is used, 
in which all dependence on renormalization or factorization scales is avoided 
by a complete resummation of ultraviolet logarithms. The Bernstein Poly- 
nomial method is used to compare these QCD predictions to the xFa data 
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of the CCFR collaboration, and direct fits for A^rJ-^ , with Nf = 5 effective 
quark fiavours , over the range 20 < < 125.9 Gev^, were performed. We 
obtain = 202^45 MeV, corresponding to the three-loop running coupling 
as{Mz) — 0.1174jlo;oo43- Including target mass corrections as well we obtain 
= 228^1 MeV, corresponding to a^iMz) = 0.11961^:^^^1. 
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1 Introduction 



The recent measurements of the CCFR collaboration provide the most pre- 
cise determination of the deep inelastic scattering (DIS) structure functions 
of neutrinos and anti- neutrinos on nucleons . In this paper we wish to com- 
pare CCFR measurements of the structure function xF-^^x^Q"^)^ with QCD 
predictions for its moments in order to determine Ajg^ g, ^, ^, Q . We intend 
to follow essentially the same method of analysis employed in Refs.]^, in 
which one constructs averages of the measured structure function with re- 
spect to suitably chosen Bernstein polynomials. The polynomials are chosen 
so that the range of x for which the experimental xF^ is not determined 
makes only a small contribution to the averages. These experimental Bern- 
stein averages are then fitted to the QCD predictions for the corresponding 
linear combinations of moments. The analysis uses the available next-next- 
to-leading order (NNLO) anomalous dimension and coefficient function re- 
sults for the odd moments with n = 1, 3, 5, . . ., 13 0. Refs.0, ^ also consider 
the F2{x,Q'^) structure function in ep DIS at NNLO. The key difference in 
our analysis is that the QCD predictions for the moments of F3 are obtained 
in the "Complete Renormalization Group Improvement" (CORGI) approach 
[0, in which all dependence on the renormalization scale , fi, and factoriza- 
tion scale , M, are eliminated. Customarily in the standard RG- improvement 
of QCD predictions for leptoproduction moments one chooses /i = xQ and 
M = yQ, with x and y undetermined dimensionless constants, so the renor- 
malization scale and factorization scale are proportional to the physical DIS 
energy scale Q. Refs. ^] make the standard choice = M = Q. In the 
CORGI improvement one instead keeps /x and M independent of Q. One 
is then forced to resum to all-orders the RG-predictable ultraviolet (UV) 
logarithms of Q , the logarithms of M and fi contained in the renormalized 
coupling constant, and the perturbative coefficients, then cancel, and one is 
left with predictions which are independent of /i and M. Crucially in this 
way one also generates the correct physical Q-dependence of the moments, 
whereas in standard improvement one omits an infinite subset of UV loga- 
rithms, so that the Q-dependence involves the unphysical parameters x and 
y. The approach is closely related to the Effective Charge formalism of Grun- 
berg [^. We find that our CORGI fits result in somewhat larger values of 
as{Mz) than those reported in Qj. 

The plan of the paper is to give a brief review of the CORGI approach 
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for leptoproduction moments in Section 2. Section 3 will contain a short 
description of the Bernstein Polynomial averages to be employed in the fits, 
and Section 4 details the results of the fits to the CORGI predictions for the 
moments. Section 5 contains a discussion and Conclusions. 



2 Leptoproduction moments in the CORGI 
approach 

Let us define the Mellin moments for the vN structure function xF3(x, Q^), 

MM)= (\^-^F,{x,Q^)dx. (1) 
Jo 

Adopting the notation of |^ we have the factorized form for the rfi^ moment, 
ca \ ^(")/^ 



Mn{Q ) = A{n) I — — ) exp(X(a)) {l+ri{n)a+r2{n)a'+r3{n)a'^+. . .) , 

(2) 

where the first three terms correspond to the operator matrix element < 
On{M) >, with M the factorization scale, and a=as{M) /t^ in terms of the 
RG-improved coupling. A{n) is an undetermined scheme-independent overall 
constant, and will be one of the parameters varied in the fits. X(a) is a 
function of the anomalous dimension coefficients di{n), 

M d < On > u ^ , / N 2 , / N S /on 

= lOn = —d{n)a — ai[n)a — a2[n)a — . . . . [6) 



<Or,> dM 



The first anomalous dimension coefficient, (i(n), |Ty] is independent of the 
factorization scheme (FS), the higher coefficients d.i{n) (i>l) define the FS. 
The coupling a satisfies the beta-function equation, 

da 

^dM ^ ^^"''^ ^ ~^^^(^ + ca + C2a^ + c^a^ + ...). (4) 

Here h = (33 — 2Nf)/6 and c = (153 — 19A^/)/126 are the first two univer- 
sal coefficients of the beta-function. The remaining coefficients are scheme- 
dependent and determine the renormalization scheme (RS). The final factor 
in Eq.(2) is the coefficient function. We use the notation d=as{f^)/7T , where 
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/i is the renormalization scale. 



The self-consistency of perturbation theory means that there is a depen- 
dence rk{n){ii, M, C2, ■ ■ ., Ck] di, d2, . . ., d^), on the parameters specifying the 
FS and RS 0. The coefficient ri depends on the factorization scale M, with, 

n(n)(M)=<i(„,)(l„f -!„£-) (5) 

Here Am„ is an FS and RS-invariant dimensionful constant associated with 
the moment, and it is the second UV logarithm in Eq.(5) which determines 
the physical Q-dependence of J^niQ^)- ^Mn is directly related to A^g^ of 
modified minimal subtraction by, 

.2c. (-c. ( di{n) ri{n)\ . . 

A«. = AWy)'-Wy;^ + ^ . (6) 



where di{n) is the MS NLO anomalous dimension coefficient, and ri(n) is 
computed in the MS scheme with M = Q. The {2c/b)~'^^'^ factor corresponds 
to the standard convention for defining A^g^ [|lT]. Using Eq.(5) one can trade 



M for ri as a parameter on which depends, similarly fx can be traded for 
fi=ri(M = /i). The self-consistency of perturbation theory allows one to 
obtain the partial derivatives of the coefficients with respect to the FS 
and RS parameters {ri, fi,C2, ■ ■ .,Ck]di,d2, . . .,dk}. On integrating these we 
obtain expressions for the as multinomials in the FS and RS parameters, 
with scheme- independent constants of integration, X^.. Thus at NNLO we 
have, 

/ \ A ^ \ 2/ N ^ / \~ / N cdi(n) d2(n) d(n)c2 „ , . 

(7) 

The constant of integration Xk{n) can only be determined given a N'^LO per- 
turbative calculation, whereas the remaining terms are "RG-predictable" and 
can be obtained from lower orders. The so-called Complete Renormalization 
Group Improvement (CORGI) result at N'^'LO corresponds to resumming to 
all-orders the RG-predictable terms, i.e. those not involving Xi{n), {i > k) 
[0]. It is easy to show that this is equivalent to working with standard RG- 
improvement in an FS and RS in which all the parameters are zero. In this 



5 



scheme one has a = a = oq, where the CORGI couphng ao(n) can be 
expressed in terms of the Lambert W function [12, |l^, 14 1, defined imphcitly 
by W{z)ex.p{Wiz)) = z, 



ao[n) 



c[l + W.^iz^iQ))] 



Here the "—1" subscript denotes the branch of the Lambert W function 



required for asymptotic freedom, the nomenclature being that of Ref.p3 



The A_A4„ is the scheme-invariant, related to by Eq.(6). One obtains the 
N'^LO CORGI result. 



d{n)/b 

1 + caQ{n) 

... + Xk{n)4{n)) . (9) 

Substituting the explicit expression for in terms of in Eq.(8) we obtain 
the NNLO CORGI result, 

Mn{Q') = A{n)[-W^,{zrXQ))r'^^^'\l + X^{n)al{n)) 

So we see that the moment is directly proportional to a power of the Lambert 
function W^i. The NNLO CORGI invariants X2{n) can be computed from 
the MS results for ri{n),r2{n),di{n),d2{n) 0. The NNLO anomalous di- 
mension coefficient d2{n) is only known for odd moments, n = 1, 3, 5, . . ., 13 
[^]. On rearranging Eq.(7) one has, 

X2(n)=\( T^yUn) ^-^^ + ^V^+ +r2(n)] . (11) 

We tabulate the resulting X2{n) CORGI invariants for n = 3, 5, . . ., 13, and 
the ratio A_M^/Ajjg determined from Eq.(6), in Table 1. We assume Nf = 5 
active quark flavours. 
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Table 1: The numerical values of the ratio 



and the CORGI invariants 



X2(n), for the odd moments n = 3, 5, . . . , 13 of xF^. 



n 


A— 


X2in) 


3 


2.268568660 


-.9283306650 


5 


2.999798808 


1.750455480 


7 


3.489368710 


3.858776378 


9 


3.870927376 


5.663749685 


11 


4.188431945 


7.264781157 


13 


4.462684796 


8.714133383 



3 Bernstein Averages of moments 

In phenomenological investigations of structure functions, for a given value 
of Q^, only a limited number of experimental points, covering a partial range 
of values of x, are available. Therefore, one cannot directly determine the 
moments. A method devised to deal with this situation is to take averages 
of the structure function weighted by suitable polynomials. We can compare 
theoretical predictions with experimental results for the Bernstein averages, 
which are defined by B 



Kk{Q^)= [ dxpnk{x)F-i{x,Q'^) 
Jo 

where Pnk{x) are modified Bernstein polynomials. 



(12) 



'T{k+\)T{n-k + l) 



X' 



2k, 



1 — X 



2\(n-k) 



(13) 



and are normalized to unity, Jq dxpnk{x) = 1. Therefore the integral (12) rep- 
resents an average of the function F^i^x) in the region Xnk — |Ax„fc<x<x„fc — 
|Ax„fc where Xnk is the average of x which is very near to the maximum of 
Pnk{x), and Axnk is the spread of Xnk- The key point is that values of F3 
outside this interval contribute little to the integral (12), as Pnk{x) decreases 
to zero very quickly. So, by suitably choosing n, k, we manage to adjust the 
region where the average is peaked to that in which we have experimental 
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data [0. Using the binomial expansion in Eq.(13), it follows that the aver- 
ages of F3 with Pnk{x) as weight functions, can be obtained in terms of odd 
moments, 

^ , (n-m(n+l) (-1)' /' <i,,«^+')«)-.)f3 (14) 

V(i + i)r(n - * + 1) £j - * - i)! io ^ ' ^ ^ 

using Eq.(l) then. 

For the NNLO QCD fits to be performed we are restricted to considering 
odd moments of xF^ for which the NNLO anomalous dimension coefficient 
d2{n) has been computed @], n = 3, 5, . . ., 13. We can only include a Bern- 
stein average, Fnk, if we have experimental points covering the whole range 
[xnk — Xnk — ^^^nk] 0, this mcaus that we can use only the 10 aver- 

ages F2(f"^)(g2), Fs^-P\Q^), fM^p\Q''), fII^p^Q^), Ft''\Q''), Ft''\Q^), 
F5^-f)(Q2), Fe^f^^HQ'), FfKQ^), and Fg^^^^HQ') i- To obtain these ex- 
perimental averages from the CCFR data for xF^ [|l|, we fit xF3(x,Q^) for 
each bin in separately, to the convenient phenomenological expression. 



X 



this form ensures zero values for XF3 at x = 0, and x = 1. A theoretical jus- 
tification of Eq.(16) may be found in Ref.||15||. Using Eq.(16) with the fitted 
values of A, B, C, one can then compute f'^^^\Q'^) using Eq.(12), in terms of 
Gamma functions. The resulting experimental Bernstein averages are plotted 
in Figure 1. The errors in the fI^^^\Q'^) correspond to allowing the CCFR 
data for xF^ to vary within the experimental error bars, including experimen- 
tal systematic errors |^. We have only included data for (5^>20GeV^, and 
our QCD fits will assume Nf = 5 active quark fiavours. This has the merit 
of simplifying the analysis by avoiding evolution through fiavour thresholds, 
whilst only reducing the number of fitted F^f^^^ points by eight. 
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Figure 1: Fit to xF^ using Bernstein averages 
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4 NNLO QCD fits to Bernstein averages for 

Using Eq.(15) the ten Bernstein averages FnkiQ"^) can be written in terms of 
odd moments M.n{Q'^), 

F2l(Q') = 7.5 (A^3(Q') - -M5(Q')) 

F3i(Q') = 13.125 (M^{Q^) - 2M^{Q^) + M^iQ^)) 
Fs2{Q') = 17.5 {M,{Q') - M^iQ')) 

F4i(Q') = 19.687 (M3(Q') - 3M5(Q') + 3M7(Q') - MglQ')) 
F42(Q') = 39.375 (A<5(Q') - 2M7(Q2) + A<9(Q')) 
F5i(g2) = 27.070 {MsiQ^) - 4M,{Q') + GMyiQ^) - Mn{Q')) 
F52(g') = 72.187 (M^iQ^) - SMriQ^) + ^Mg{Q^) - 4Mg{Q^)) 
F6i(g') = 35.191 (M3(Q') - 5M5(Q') + 10M7(Q') - M9(Q') 

+ 5Mii(Q')-A4i3(Q')) 
F62(g') = 117.30 (M^iQ^) - iMriQ^) + QM^iQ^) - iMuiQ^) 

FmiQ') = 187.69 {MriQ^) - 3Mg{Q') + 3Mn{Q') - Mis{Q')),17) 

We shall use the NNLO CORGI result of Eq.(lO) for the QCD prediction 
of A4n{Q'^)- The basic fit parameters will be the unknown normalization 
constants A{n), n = 3, 5, 7, . . . , 13, and A-^ , related to the CORGI A_a^„ 
by Eq.(6), see Table 1. Thus there are 7 parameters to be simultaneously 
fitted to the experimental Fnk{Q^) averages. Defining a global foi" 
the experimental data points of Figure 1 , we found an acceptable fit with 
minimum x^/d.o.f. = 1.2905/43. The best fit is indicated by the curves in 
Figure 1. Allowing within 1 of the minimum to estimate an error gives, 

Ag^ = 2021^^ MeV , (18) 

which corresponds to the three-loop MS coupling at the Z-mass, 

a,{Mz) = 0.1174l°:°°ll . (19) 
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The minimum values for all 7 fitting parameters are tabulated in Table 2. 
To attempt to include target mass corrections (TMC) in the fits we amended 
the expression for J^niQ"^) to |TH], 



2(. + 2) g^^"-^^^) + ^(#- ^20) 

2 

This results from the series expansion in powers of of the Nachtmann 
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r m„ 



moments IT^]. Here nip is the proton mass. The influence of the 0{-^) terms 
is very small, and we shall neglect them. Minimizing with the amended 
expression for the moments then gives a best fit of comparable quality, with 

^TMC) ^ 228+1 MeV . (21) 
Corresponding to the coupling at the Z-mass, 

a.(Mz) = 0.1196l°;°03l- (22) 

The best fit values of the 7 fitting parameters including TMC, are tabulated in 
Table 3. We can compare this value of as{Mz) with the corrresponding values 
as{Mz) = 0.1153±0.0041 obtained in Ref.||, and ^^(Mz) = 0.1187±0.0026 
obtained in Ref . [Q . Finally we can confirm the expectation that at the energy 
scales Q2>20GeV^ included in our fit, higher-twist (HT) effects should be 
small. We modified the expression for A4n{Q^) by an additional term p|, 

-MrHQ') = ^ (^) -Mn(Q^) , (23) 

where p is an additional phenomenological parameter which will be fitted to 
the data. The best 8 fit parameters are tabulated in Table 4. The inclusion 
of HT terms shifts the central value of by only 2 Mev (cf. Table 2) , 
confirming as expected that HT effects are negligible. Our fitted value of the 
HT parameter p = — 0.77±0.23 is to be compared with the value — 0.14±0.6 
obtained in Ref.0, and -0.31±0.80 obtained in Ref.|[. 

We should stress that the errors in the values of as{Mz) quoted in 
Eqs.(19),(22), reflect the errors in the F^^^^'^ values in Fig. 1 , to which 
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the NNLO CORGI predictions for the moments of Eq.(lO) have been fit- 
ted. In the CORGI approach all dependence on the unphysical factorization 
and renormalization scales is eliminated by the complete resummation of 
RG-predictable UV logarithms, the remaining uncertainty in the QCD pre- 
diction then resides in the unknown N^LO CORGI invariant X^ln). It would 
in principle be straightforward to use Pade approximants to estimate the un- 
known d^ln) anomalous dimension coefficients, as in Ref.[Q, and perform fits 
with an estimated X^ln), but we have not done so in this work. 

Table 2:Numerical values of fitting parameters, for the best fit of Figure 1. 





MeV 


A{3) 




A{5) 


n 1 c;Q-tJ-uo8 
u-J- jy+0.008 


A{7) 


n n7-u.uui3 

^•^ ' +0.0012 


A{9) 


n n.+O.UOU 
U-U4_0.0013 


Aill) 


n nQ+u.ooio 


A(13) 


n, ,^90-0.0020 
U.UZJ^Q 0019 



Table 3:Numerical values of fitting parameters, including TMC. 





228+;^^ MeV 


A{3) 


0.48lJ:Ui^ 


A{5) 


U--LO+0.004 


A{7) 


n ns-U 0024 

U-UO^0.0021 


A{9) 


U-UO+0.0027 


A{11) 


n nq-u.uoas 

U-UO+0.0033 


A(13) 


u.uuy+0.0029 
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Table 4:Numerical values of fitting parameters, including HT. 





204:+^^ MeV 


A{3) 




A(5) 


U-J-Oy+0.008 


A{7) 


n n7-U-00i3 

' +0.0012 


A{9) 


Q Q4+0.U0i4 
U-U4_0.0013 


A{11) 


n nQ+u.oolo 


A(13) 


U.UZy^Q QQ29 


P 


n 77+0.^3 
— U. 1 ^_o.23 



5 Discussion and Conclusions 

In this paper we liave used a similar method of analysis to that of P] to fit 
QCD predictions for the moments of the uN DIS structure function xF^, to 
suitably constructed Bernstein polynomial averages of the CCFR experimen- 
tal data for xF^ The key difference in our approach has been the use 
of Complete Renormalization Group Improved (CORGI) 0, | NNLO QCD 
predictions, hence avoiding the need to make particular ad hoc choices of the 
dimensionful renormalization scale, /i, and factorization scale, M. The most 
important motivation for the CORGI approach is that by completely resum- 
ming all the UV logarithms one correctly generates the physical dependence 
of the moments A^n(Q^) on the DIS energy scale Q. From Eq.(lO) this is 
seen to be. 



M„(g2)~A(n) 




MnJ 



-d(n)/b 



b\n{Q/A 



Mr 



(24) 

so that the large-Q behaviour is controlled by the Lambert W^i function, 
and the ratio Q/A_x4^, with A_m^ the FS and RS- invariant defined in Eq.(5). 
In contrast with standard RG-improvement and choosing /i = M = xQ, with 
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X a dimensionless constant, one has instead, 




) 



b/c\ 1 ~d(n)/h 



( 



( 



1 



» 

(25) 



-1 



1 + 



b\n{xQ/Aj^) 



which manifestly depends on the unphysical parameter x. Here Aj^={2c/b) ^ 
xAjjg. In Refs.p, |^, |] the standard choice a; = 1 is made. 

It seems clear to us that CORGI QCD predictions should be used in 
the fits. However, the value of as{Mz) = 0.1196±0.003 that we obtain, in- 
cluding TMC effects, is consistent with the result as{Mz) = 0.1153±0.0063 
obtained in Ref. ||^ using a similar method of analysis, but with the stan- 
dard X = 1 scale choices. We are in even closer agreement with the value 
as{Mz) = 0.1195±0.004 reported in Ref.[§|, but these authors use a very 
different method of analysis involving the Jacobi Polynomial technique, to- 
gether with the standard x = 1 scale choice. This indicates that the incom- 
plete resummation of UV logarithms implicit in the standard scale choice 
IJ, = M = Q, does not greatly modify the fit. This is underwritten by the 
reasonably small NNLO CORGI invariants X2{n), and ratios A^^/A^g^, ap- 
pearing in Table 1. From Eq.(6) and Table 1 one can see that the CORGI re- 
sult corresponds to using an n-dependent scale choice, ranging from x = 0.44 
for the n = 3 moment to a; = 0.22 for the n = 13 moment. The CORGI 
approach also differs in that C2 = and di{n) = d2{n) = 0, rather than the 
MS values. It should be noted, however, that use of the standard physical 
choice of renormalization scale is in general likely to result in misleading de- 
terminations of A;^, as has been discussed in Ref.[§] in the case of e~^e~ jet 
observables. 
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